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We study analytically and numerically the winding of a flux line around a columnar defect. 
Reflecting and absorbing boundary conditions apply to marginal or repulsive defects, respectively. 
In both cases, the winding angle distribution decays exponentially for large angles, with a decay 
constant depending only on the boundary condition, but not on microscopic features. New non- 
universal distributions are encountered for chiral defects which preferentially twist the flux line in one 
direction. The resulting asymmetric distributions have decay constants that depend on the degree of 
chirality. In particular, strong chirality encourages entanglements and leads to broad distributions. 
We also examine the windings of flux lines in the presence of point impurities (random bonds). Our 
results suggest that pinning to impurities reduces entanglements, leading to a narrow (Gaussian) 
distribution. 
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I. INTRODUCTION AND SUMMARY 



Winding angles of paths are of great interest not only 
in mathematics, but also in the physics of polymers, and 
flux lines in high-T c superconductors. The topological 
constraints produced by the windings of polymers B or 
magnetic flux lines ]||| around each other, result in en- 
tangled phases with slow dynamics. The simplest case 
that can be studied is the winding of a 2-dimensional 
random walk around a point, or equivalently, a flux line 
in 3 dimensions around a columnar pin H]. In 1958, 
Spitzer [g| showed that the probability distribution for 
the winding angle 9 of a Brownian path around a point 
is a Cauchy law for large time i, i.e. 
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Similar Cauchy laws are obtained for winding around 
several points in 2 dimensions and around several 
straight lines in 3 dimensions 0. These results are ob- 
tained by employing a variety of techniques such as stan- 
dard diffusion equation J |,||J^] , path integrals |]|l(]], or 
probability theory pl|,|s|j7j. By contrast, the winding an- 
gle of a self-avoiding walk in 2 dimensions obeys a Gaus- 
sian distribution, the scaling variable being x — #/V '4 lnt 
(See also Ref. [|l3| for an expansion around four di- 
mensions.) 

As pointed out in Refs. jl4|,[l5|, the above Cauchy law 
has pathological properties which make its relevance to 
any physical situation questionable. In particular, be- 
cause of the slowly decaying tails at large x, the averages 
of both 8 2 and \9\ are infinite. The origin of this diver- 
gence is that a finite segment of the Brownian walk can 
wind infinitely often around a point center. While this 
is correct for an idealized random walk, in any physical 
system one expects a cutoff due to either finite diameters 
or stiffness. The case of a Brownian walk in 2 dimensions 
around a disc of finite diameter was studied in Ref. . 
The resulting winding angle distribution is |l4| , |il| 
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In Ref. fig] , the following result for the winding angle 
distribution for a random walk with steps of finite size is 
derived 
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The same result is obtained in Ref. || for the distribution 
of "big windings" of Brownian motion around two point- 
like winding centers. Saleur Jlq | suggests that the differ- 
ence between Eqs. (d) and @is due to different bound- 
ary conditions for the walkers at the winding center. A 
review of many topological and entanglement properties 
of polymers can be found in Ref. fDj ], 

In this paper, we further study the issue of winding 
angle distributions and their universality, with particu- 
lar emphasis on their applicability to magnetic flux lines 
(FLs) in high-Tc superconductors. We start in Sec. [n]by 
providing a derivation of Eqs. (|l|)-(||) based on conformal 
properties of random walks that does not require any ad- 
vanced mathematical techniques. It illustrates well the 
origin and universality of the exponential tail in Eqs. (|) 
and (jj^), and explains the factor 2 between their two de- 
cay constants. We argue that these two cases are appli- 
cable respectively to the windings of a flux line around 
a repulsive or marginal columnar defect. Actually, most 
columnar defects are attractive and localize the flux line 
to their vicinity. The corresponding probability distribu- 
tions (for walkers with initial and final poin ts clo se to the 



winding center) are also calculated in Sec. [I D 



To test the generality of the analytical results, we per- 
formed a number of numerical tests that are described 



in Sec. [II. Simulations of random walks were performed 



on both square and cubic lattices. For reflecting bound- 
ary conditions we chose a winding center that was shifted 
from the lattice sites crossed by the walker. For absorb- 
ing boundary conditions, the center was one of the lat- 
tice sites which the walkers were not allowed to cross. A 
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transfer matrix method can be used to evolve the wind- 
ing angle distribution with the number of steps t. De- 
spite a rather slow convergence to the asymptotic limit, 
the numerical results do indeed support the universality 
of distributions given in Eqs. (^) and (||). 

In the course of numerical tests, we encountered one 
case (reflecting boundary conditions for a directed path 
along the diagonal of a cubic lattice) which did not con- 
form to any of the above expected universality classes. 
Further examination revealed that we had inadvertently 
constructed a chiral defect: upon encountering the defect, 
the path had a statistical advantage to wind in one direc- 
tion as opposed to the other. In the language of random 
walks, this corresponds to a rotating winding center. The 
breaking of the symmetry at the center is in fact a rele- 
vant perturbation, and leads to a new type of probability 
distribution as discussed in Sec. 0. Although we do not 
yet have a complete analytic understanding of this novel 
universality class, we can account for certain features of 
this distribution. In particular, weak chirality leads to 
narrow asymmetric distributions, while for strong chiral- 
ity both wings of the distribution are widened, i.e. strong 
chirality of the defect enhances entanglements. 

All the distributions examined in Sees. pl^Ty] describe 
the windings of ideal random walks, and in all cases 
the appropriate scaling variable is the combination x = 
28 /hit. This universal feature is due to the Markovian 
nature of these walks, and can be explained as follows: 
After time t, the walker has a typical distance r(t) cx \fi 
from the starting point, which is chosen to be close to the 
winding center. Assuming that r(t) is the only relevant 
length scale, dimensional arguments, combined with the 
Markovian property, suggest that dr/dO = rf(9). The ro- 
tational invariance of the system implies that f(9) must 
be a constant, i.e. the increase in winding angle cannot 
depend on the number of windings or angular position. 
Hence, 



dr dt 
do oc — cx — = fl(lnt) 
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leading to a scaling variable proportional to 9 /hit. 

Because of their non-Markovian nature, we cannot ap- 
ply the arguments of the preceding paragraph to self- 
avoiding walks. Indeed, such walks have a Gaussian 
winding angle distribution in the scaling variable x = 
0/y/hit |Q. The self-similarity of the walk suggests the 
following scaling argument: Starting from the origin di- 
vide the walk into segments of 1, 2, • • • , 2™ w t/2 steps. 
Since the a th segment is at a distance of roughly 2 av from 
the center (with v — 3/4) and has a characteristic size 
of the same order, it is reasonable to assume that each 
segment spans a random angle 9 a of order one. Under 
the mild assumption that the sum 9 — J2a=i ®a satis- 
fies the central limit theorem, we then conclude that 9 
is Gaussian distributed with a variance proportional to 
n oc In t. 

We can ask why the above argument of self-similarity 
does not apply to ideal walks. The reason is the relevance 



of the finite winding center. An important difference be- 
tween ideal and self-avoiding walks is that the probability 
of returning to the winding center at the origin asymp- 
totically vanishes for the latter because of the exclusion 
zone set up by self-avoidance. It is thus expected that 
properties of the winding center (size, absorbing versus 
reflecting nature, chirality) are irrelevant for self- avoiding 
walks. The ideal walks, on the other hand return to the 
origin quite often, and upon rescaling see winding cen- 
ters of different size. The assumption that different scaled 
portions of the ideal walk are self-similar is not correct. 

There is one case where both arguments should hold: 
An ideal walk around a point winding center. The argu- 
ment based on self-similarity actually states that the final 
distribution is obtained from the composition of In t in- 
dependent random variables. If each variable has a finite 
variance, the overall distribution will be Gaussian. If not, 
other (Levy) distributions are possible. The Cauchy dis- 
tribution is in fact a limiting case for widely distributed 
variables. The requirement that both arguments should 
hold, immediately selects a Cauchy distribution! 

The scaling argument, which by no means is claimed to 
be exact, should apply to other self-similar walks where 
the probability of return to the origin is small. An inter- 
esting example is provided by directed paths in random 
media Jl8|], which for example describe the behavior of 
a flux line in the presence of (quenched) point impuri- 
ties. Typical wanderings of such paths scale as t", with 
v « 0.59 > 1/2. The pinning by impurities greatly re- 
duces the probability of the walker returning to the ori- 
gin, and the above arguments again suggests that the 
winding angle distribution is Gaussian. Although cer- 
tainly not conclusive, the numerical results reported in 
Sec. |v| appear to support this prediction. Typical wind- 
ing angles thus scale as Vint as opposed to lni in the 
pure case. This simple example thus suggests that topo- 
logical entanglements become relatively less important in 
the presence of pinning to impurities. 

The results of this paper point out the rich behav- 
ior already present in the simplest of problems involving 
topological defects. Properties of the winding center (fi- 
nite size, chirality, • • •), interactions, various types of ran- 
domness are all potentially relevant, leading to different 
universal distribution functions. These results could also 
produce some interesting physical manifestations. For 
example, we demonstrate in Sec. VI that there is a sharp 



crossover between free and coiled configurations, if there 
is an energy proportional to the winding number. 



II. DERIVATION OF WINDING ANGLE 
DISTRIBUTIONS 

A. Spitzer's law 

Before studying the winding of a Brownian path 
around a center of finite diameter, we first sketch a 
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derivation of Spitzer's law in Eq. (|y). We follow the ap- 
proach in Ref. jllj, translating it into a more physical 
language that does not rely on a familiarity with martin- 
gales. A basic ingredient is the invariance of Brownian 
motion under conformal mappings. Let 



z(t) = xi(t) + ix2(t) 



(5) 



represent the original walk in the complexified two- 
dimensional plane. The time evolution of each random 
walker satisfies 



dz = r)(t)dt, 
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where the complex random velocity has zero mean and 
is uncorrelated at different times, with 



{r,{t)r)*{t'))=2D5{t-t'). 
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The radius of the walker, and its winding angle can be 
extracted from 



C(t) = \nz{t) = p{t)+i6(t), 



(8) 



where p = In r = In \fx\ 



the stochastic motion of the walker in the new complex 
plane is highly correlated to its location. This feature 
can be removed by defining a new time variable 



dr = 



dt 



\z(t)\ 2 

for each walker, which leads to 



dC, = p(r)dr, with p(t) = z*(t)i](t). 



(9) 
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Since 



(p(t)p*(t')) = 2D\z(t)\ 2 S(t - t') = 2DS(t - r'), (11) 

the evolution of £(t) is that of a Brownian walk. 

For simplicity we choose the initial condition £(t = r = 
0) = 0, i.e. the original walker starts out at z = 1. We 
also set the diffusion constant to D = 1/2, so that the 
mean square distance over which the walk moves during 
a time t is (r 2 (t)) = t. Since there is a separate transfor- 
mation r(t) for each walker, it is not useful to subdivide 
the walkers £ according to r. Instead the walkers of in- 
terest are the ones that reach lines of fixed p for the first 
time, independent of the value of r. This is because af- 
ter a time t, the value of p is almost certainly equal to 
Int. In fact, the probability that r(t) is within an interval 
[v 7 ^ 1-6 ^ 2 , \fTTt^ 1+e ^ 2 ] around its mean value of ypUt, 



p(t,e) 



v^t (1+t)/2 



exp (-r 2 /2t) 
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approaches unity in the limit t — > oo. In this limit, 
the distance r from the starting point z = 1 is iden- 
tical to the distance from the origin, and p(t, e) is 
identical to the probability that £(r) is in the interval 
[0.5(1 -e) hit, 0.5(1 + e) hit]. 

We next shrink the complex plane £ by a factor of 
(lni)/2. The walkers of interest are now those that hit 
the line with real value of unity for the first time. The 
imaginary coordinate x, of the hit is now related to the 
winding angle by x = 29(t)/lnt. Thus, in the limit 
t — ► oo, the probability p(x)dx is precisely the proba- 
bility that a Brownian walker (2£/ \nt) starting from the 
origin, hits a vertical line at a distance one from the ori- 
gin, for the first time, at a height between x and x + dx. 

Since the walkers hit the line for the first time at dif- 
ferent times r, we need the probability density p\(r) that 
the first hit is at t. Its integral -Pl(t) = f Q dr'p^T 1 ), is 
the probability for the walker to hit the line at least once 
before time r. Since after hitting the wall the walker is 
equally likely to proceed in either direction, the latter is 
also twice the probability that the endpoint of the walk 
is beyond the line at r, and hence given by 



Since d( = 7](t)dt/z(t), Pi(r) = 2 
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The derivative gives 
dP x 

Pl{T) = ^ 



exp (-r 2 /2r)_ dr = r exp(-l/2s) ^ 
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Since the vertical and horizontal components of the walk 
are independent, we finally obtain 



p(x) 



00 , exp (-l/2r) exp (-x 2 /2t) 
/ — ' 
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which is Spitzer's law; exact in the limit t — > oo. 



B. Repulsive defects and finite absorbing centers 

The unphysical situation that a path element of finite 
length can wind infinitely often around the origin can be 
removed by introducing a winding center of finite size. 
A FL winding around a repulsive columnar defect can- 
not penetrate the defect, and all configurations where 
the FL and the defect intersect have to be removed [jlJJ . 
This situation corresponds to a random walk around an 
absorbing winding center. Two directed polymers with 
hard-core interaction winding around each other are de- 
scribed by the same situation, the radius r now being the 
relative distance between the two polymers. 

We choose a disc of radius R < 1 around the origin 
which cannot be entered by the Brownian walk z(t). A 
path of length t now has a maximum winding angle of 
t/R. The walk is again started from z(0) = 1, and after a 
long time t, its end point is almost certainly at a distance 
r between y / 7r^ 1_e -'/ 2 and y/nt^ 1+e ^ 2 from the origin, as 
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given by Eq. ( |l2| ) . (The excluded disc of radius R <C \fi 
has negligble influence in this limit.) The random walks 
in the complex plane C = l nz consequently have real 
values in the interval [0.5(1 — e) ln£, 0.5(1 + e) ln£] as be- 
fore. But now the walks £(t) cannot go to the left of 
a line at a = \nR. The probability density pa{x) for 
having a scaled winding angle x = 29/\nt is therefore 
identical to the probability density that a random walk 
starting at the origin hits a line at distance one, for the 
first time, at height x without going beyond a line at 
distance 2| lni?|/lnf on the opposite side. 

Let P 0j b(y,r) denote the probability that a Brownian 
walk starting at y £ [a,b] hits the point b for the first time 
during a time interval r, without hitting point a previ- 
ously. If we consider the points a and b as absorbing 
boundaries, this is identical to the probability that the 
walk is absorbed at boundary b before time r. Since for 
sufficiently small At, the walker is only a short distance 
Ay from its starting point, we have 



Pa,b{y,r) = f 
Jo 

X [Pa.biV 



d(Ay) 
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Ay, 
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Expanding the above equation to the order of At indi- 
cates that P a ,b{y,T) satisfies a diffusion equation. The 
appropriate boundary conditions are P a j,{a,T) = and 
Pa,b(b, t) = 1 with the initial value P a ,b(,y, 0) = 0, result- 
ing in M 



Pa,b(y,r) = 



y- a 
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x exp 
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The probability that the walk is absorbed at the right- 
hand boundary during the time interval [r, r 4- dr] , is 
dTd T P at b(y,T). Note, however, that / °° dTd T P a: b(y, t) — 
(y — a) I (b — a) , i.e. equal to the total fraction of particles 
absorbed at the right-hand boundary (inversely propor- 
tional to the separations from the boundaries). To calcu- 
late pa(x), we need the fraction of these walks absorbed 
between r and r + dr, equal to ((6 — a)/(y — a))d T P. 
Hence (with a = 2 lni?/ lni, b — 1 and y = 0) 
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The last step is achieved by first performing a Fourier 
transform with respect to x, followed by integrating over 



r, and finally inverting the Fourier transform. (Alterna- 
tively the t integration can be performed by the saddle 
point method.) In the limit of large t, the variable a is 
very small, and we can replace the sine-function by its 
argument. Taking the sum over v, we find 



Pa(x) = 



exp [ttx/ (1 — a)] 



l-a{exp [ttx/(1 - a)} + l} 2 



(13) 



Changing the variable from x to 
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' (l-o) ln(i/i? 2 )' 

and noting that pA(x)dx = pA(x)dx, leads from Eq. ( p^ ) 
to 



PA x 



26* 



ln(t/R 2 



4 cosh (ttx/2) 
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This is identical to Eq. (||), since x = x for large i. The 
use of the variable x is, however, more appropriate, since 
it makes the argument of the logarithm dimensionless, 
and explicitly indicates the basic time unit. The above 
distribution, which is exact in the limit t — > oo, has an 
exponential decay for large x, in contrast to Spitzer's law 
in Eq. ([!]). In particular, the mean winding angle is now 
finite. The analogy to random walkers in the plane £, 
confined by the two walls, provides simple physical jus- 
tifications for the behavior of the winding angle. In the 
presence of both walls, the diffusing particle is confined 
to a strip, and loses any memory of its starting position 
at long times. The probability that a particle that has 
already traveled a distance 9 in the vertical direction pro- 
ceeds a further distance d9 without hitting either wall is 
thus independent of 9, leading to the exponential decay. 
On the other hand, if there is no confining wall on the 
left-hand side, the particles may diffuse arbitrarily far in 
that direction, making it less probable to hit the wall on 
the right hand side. The absence of a characteristic con- 
fining length thus leads to the unphysical divergence in 
Spitzer's law. 

We also expect that different shapes of the winding 
center do not modify Eq. (p^). This is because of the 
division by ln£/2, which maps any excluded shape to the 
vertical line at a — for long enough times. However, 
different shapes should result in different effective values 
for R, and the time constant in Eq. (|l4|). For long times, 
the polar angle of the particle may take any value be- 
tween and 2ir with equal probability, and the particle 
sees an averaged radius of the winding center. 

The crossover to Spitzer's law in the limit R — > is not 
apparent from the above expression for pa{x). The tran- 
sition can be made in the initial equation for pa{x), but 
not in the final result of Eq. (|lj) , where time is measured 
in units of R 2 . The time scale in Spitzer's law is set by 
the distance of the initial point from the winding center, 
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which is infinite in units of the radius of the winding cen- 
ter. Consequently one time unit in Eq. (jl|) corresponds 
to infinitely many time units in E q. ([Til ), and the limits 
i^ooori?^0ori-^0in Eq. ( |14f ) correspond to the 
limit t — ► or x — ► oo in Eq. ([!]) . 



C. Neutral defects and reflecting winding centers 



For windings around several point centers, or the wind- 
ing of a random walk on a lattice around a point different 
from the vertices of the lattice pjlq], no configurations 
are forbidden, and no walks are removed. The result- 
ing winding angle distribution in both cases is given by 
Eq. The random walk on a lattice can be regarded 
as a model for a directed polymer of finite stiffness, the 
step size being of the order of the persistence length. 
This situation applies to magnetic FLs winding around a 
columnar defect that is marginal, i.e. neither attractive 
nor repulsive. We can obtain Eq. ([}]) by repeating the 
calculations of the previous subsection, but replacing the 
absorbing boundary condition P a ^(a, r) = with the re- 
flecting condition dP ai b(y,T)/dy \ y=a = 0. There is thus 
no current leaving the system at point a, and walkers 
which hit the winding center are reflected. We then find 



P a ,b{y,r) = 1- 
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y — 



/2 
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ir{v+\/2){b-y) 



x exp 
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2 



and 



Pr(x) 



1 fir (v + 1/2) 
b — a 



dP o ,i(0,r)exp (-x 2 /2t) 



(15) 




u=0 



x exp 
1 



a \ 1 — a 

tt(i/ + 1/2)|x| 



{I -a) 



2 cosh(7ra/2) 



(16) 



where again time is measured in units of R 2 , and the limit 
t — > oo has been taken. For large x, where the walk has 
lost the memory of its initial distance from both walls, 
this probability decays exponentially as exp [— %x/2]. A 
random walk confined between an absorbing and a re- 
flecting wall that have a distance 1 can be mapped to 
a random walk confined between two absorbing walls at 
distance 2. After rescaling the wall distance and the x- 
coordinate by two, this explains the factor 1/2 between 
the decay constants in the tails of the distributions in 
Eqs. (|) and (§). 



D. Attractive defects and returning walks 



A magnetic FL winding around an attractive colum- 
nar defect of radius bo and binding energy Ug per unit 
length, is bound to that defect. If the temperature is 
above a crossover temperature T* cx bo^/TJo, the line is 
only weakly bound and wanders horizontally over the dis- 
tance of the localization length lj_(t) ~ bo exp[(T/T*) 2 ] 
[Q. The mean vertical distance l z between consecutive 
intersections of the FL with the defect is consequently 
proportional to l\. Over this distance the FL can be 
approximated by a directed walk which returns to its 
starting point (the winding center) after a time l z . 

Using the result in Eq. ([bfl), we can derive the wind- 
ing angle distribution p\{x) for such confined random 
walks. Each walk that returns to its starting point af- 
ter time t (in the z-plane) is composed of two walks of 
length t/2 going from th e sta rting point to z(t/2). As we 
have seen in subsection [II A| , almost all walks of length 
t/2, when mapped on the plane 2£/ln(t/2), have their 
endpoint on a vertical wall at distance 1 from the origin. 
The winding angle distribution for these walks is given by 
Eq. (Ibj), with t replaced by t/2. The probability that a 
walk that returns to its starting point has a winding angle 
9 is therefore obtained by adding the probabilities of all 
combinations of two walks of length t/2 whose winding 
angles add up to 8, i.e. 



p° A (x) = dy 



4 cosh (it y) A cosh (tt(x — y)) 



, „ coth(— )- 1 , (17) 
2sinh 2 (^i/2) V 2 2 / 

where x = 29/ \n(t/2R 2 ). Surprisingly, this result is dif- 
ferent from the one given in Eq. (4.1) of Ref. |L6]]. We are 
not sure of the origin of this discrepancy, but note that 
the result in Ref. U& has a divergence at x = which is 
unphysical. 

If in place of convoluting two absorbing probability 
distributions, as in Eq. (|l7|), we use the probability dis- 
tributions for reflecting boundary conditions, the final 
result is modified to 



2 sinh(7nr/2) 



(18) 



For large x w x, the above two expressions decay as 
xexp[— 7rx] and a; exp [— ttx/2] respectively. A FL of 
length T is roughly broken up into T/l z segments be- 
tween contacts with the attractive columnar defect. We 
can assume that the winding angle of each segment is 
independently taken from the probability distribution in 
Eq. (17) with t sa l z . Adding the winding angle dis- 
tributions of all segments leads a Gaussian distribution 
centered around 5 = 0, and with a variance proportional 
to Tln(l z )/l z . 
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III. UNIVERSALITY OF THE WINDING ANGLE 
DISTRIBUTION 



At first sight it may seem surprising that the winding 
angle distribution in Eq. ( |l6| ) for random walks around 
a reflecting center of finite diameter is identical to the 
distribution of "large windings" for Brownian motion 
around two point-like centers ||, or to distribution for 
random walks of finite step size p5[ . The increase in 
winding angle is largest when the walk is close to the 
winding center, and should thus be quite sensitive to the 
lattice structure, or the ratio between step size and wind- 
ing center size. However, a closer examination of the 
random walks in the scaled £-plane reveals that, in the 
limit t — > oo, almost no increase in the variable x occurs 
when the walker is within a certain finite distance from 
the winding center (corresponding to an infinitely small 
distance from the reflecting wall in the scaled C-plane). 
Hence, microscopic details do not enter the winding an- 
gle distribution, and Eq. (16) holds as long as there is a 
winding center of nonvanishing size (or, equivalently, an 
upper cutoff to the maximum possible winding angle per 
unit time), and as long as no walk coming close to the 
winding center is absorbed. 

If a walk hitting the winding center is absorbed with a 
nonvanishing probability, almost all walks are absorbed 
in the limit t — > oo, and the winding angle distribution of 
the surviving walks is given by Eq. ( I4j) . As for reflecting 
boundary conditions, microscopic details of the system 
do not influence this distribution. Apart from these two 
universality classes, we shall later encounter a third class 
of walks corresponding to random motion around a ro- 
tating winding center; the winding angle distribution in 
this case depends on the rotation speed of the center. In 
this section, we check numerically the former situations, 
where the symmetry between ±9 is not broken. 

We performed computer simulations of random walks 
on a square lattice with both reflecting and absorbing 
boundary conditions. Reflecting boundary conditions are 
simply realized by choosing a winding center different 
from the vertices of the lattice, and thus never crossed 
by the walker. On the other hand, to simulate absorbing 
boundary conditions, the winding center is chosen as one 
of the lattice sites (say the origin), but no walk is allowed 
to go through this point. 

The winding angle distributions are most readily ob- 
tained using a transfer matrix method which calculates 
the number of all walks with given winding angle and 
given endpoint after t steps, from the same information 
after t — 1 steps. The winding center is at (0.5, 0.5) for 
reflecting boundary conditions, and at the origin for ab- 
sorbing boundary conditions. The walker starts at (1,0), 
and the winding angle is increased or decreased by 2tt 
every time it crosses the positive branch of the Xi-axis. 
Due to limitations in computer memory, we applied a cut- 
off in system size and winding angle for times t > 120, 
making sure that the results are not affected by this ap- 



proximation. The largest times used, t = 9728, required 
approximately 3 days to run on a Silicon Graphics Indy 
Workstation. 
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FIG. 1. Winding angle distribution for random walks on a 
square lattice with reflecting boundary conditions for t — 38 
(dot-dashed), 152 (long dashed), 608 (dashed), 2432 (dotted), 
and 9728 (solid). The horizontal axis is x = 26»/ln(2f). The 
thick solid line is the analytical result of Eq. (0) . 



Figures [l] and |^ show the simulation results for the 
two cases. The asymptotic exponential tails predicted 
by theory can clearly be seen; deviations from the the- 
oretical curve for smaller values of the scaling variable 
x = 29/ ln(2i) are due to the slow convergence to the 
asymptotic limit. Since the scaling variable depends log- 
arithmically on time, the asymptotic limit is reached only 
for large hit. Note that the only free parameter in fitting 
to the analytical form is the characteristic time scale ap- 
pearing inside the logarithm. With t measured in units of 
single steps on the lattice, we found that a factor 2 in the 
scaling variable provides the best fit. In the limit t — > oo, 
different scales of t give of course the same asymptotic 
winding angle distribution. 
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FIG. 2. Winding angle distribution for random walks on a 
square lattice with absorbing boundary conditions. The sym- 
bols and the variable x are the same as in the previous figure. 
The thick solid line is the analytical result of Eq. (^) . 
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To further test the universality of these distributions, 
we also performed simulations of a flux line (directed 
path) proceeding along the diagonal of a cubic lattice 
in three dimensions. The FL starts at (1,0,0), and at 
each step increases one of its three coordinates by 1. We 
determined the winding angle distribution around the di- 
agonal (1,1,1)- direction, excluding from the walk all 
points which are on this diagonal (a repulsive columnar 
defect, corresponding to the case of absorbing winding 
center). The excluded points lie on the origin when the 
FL is projected in a plane perpendicular the diagonal. A 
cutoff of 243 in system size was imposed for the trans- 
fer matrix calculations. The winding angle distribution 
p(x = 29/ hut) is shown in Fig. || for different times. As 
for the square lattice, an exponential tail with decay con- 
stant of 7T can be clearly seen. Our numerical results, as 
well as the analytical considerations, thus indicate clearly 
that the winding angle distributions for reflecting and 
absorbing boundary conditions are universal and do not 
depend on microscopic details. 

Due to the special properties of directed paths along 
the diagonal of the cube, the case of reflecting boundary 
conditions leads to an asymmetry between windings in 
positive and negative directions. This is because it takes 
only three steps to make the smallest possible winding 
in one direction, but six steps in the opposite direction. 
This situation is discussed in detail in Sec. 0. 




30.0 



FIG. 3. Winding angle distribution around the preferred 
direction for a flux line (directed path) in 3 dimensions for 
t = 243 (dot-dashed), 729 (long dashed), 2187 (dashed), 
6561 (dotted), and 19684 (solid). The horizontal axis is 
x = 2(9/ln(2i). The thick solid line is the analytical result 
ofEq. (0). 



IV. CHIRAL DEFECTS AND ROTATING 
WINDING CENTERS 



So far, we only considered situations that are symmet- 
ric with respect to the angles ±0. For directed paths 
on certain lattices, however, this symmetry is broken, as 
mentioned in Sec. Ill, A directed walk that proceeds at 
each step along +xi, +X2, or +X3 direction on a cubic 



lattice can be mapped on a random walker on a two- 
dimensional triangular lattice as indicated in figure (|(a). 
Each bond can be crossed only in one direction, and the 
winding center for reflecting boundary conditions must 
be different from the vertices of the lattice, ft is apparent 
from this figure that the random walker can go around 
the center in 3 steps in one angular direction, but in no 
less than 6 steps in the other direction. 




a) 



b) 



FIG. 4. Triangular and square lattices with directed bonds. 
The winding centers are indicated by a circle o. 

An alternative description is obtained by examining 
the position of the walker after every three time steps. 
The resulting coarse-grained random walk takes place on 
a regular triangular lattice, but now the walker has a fi- 
nite probability of 3 x 2/3 3 = 2/9 of staying at the same 
site. If this site is one of the three points next to the 
winding center, the winding angle is increased by 27r in 
one of the six possible configurations that return to the 
site after three steps. In other words, the walker has a 
finite probability of having its winding angle increased in 
the proximity of the center. The amount of this biased 
increase in angle depends on the structure of the lat- 
tice and will be different for other directed lattices. An 
equivalent physical situation occurs for Brownian motion 
around a rotating winding center, e.g. a rotating reflect- 
ing disc which does not set the surrounding gas or liquid 
into motion. 

Since this angular symmetry breaking is already 
present in the above simple example of a directed walk, it 
is also likely to occur in more realistic physical systems, 
such as with screw dislocations. We thus use the term 
chiral columnar defect, to indicate that each time the FL 
comes close to the defect, it finds it easier to wind around 
in one direction as opposed to the other. (Of course, to 
respect the reflecting boundary conditions, there must be 
no additional interaction with the defect.) 

After mapping to the rescaled C-plane introduced in 
Sec. [n| the above situations can be modeled by an up- 
ward moving reflecting wall on the vertical axis. Each 
time a Brownian walker hits this wall, its vertical posi- 
tion x = 20/ In i is increased by a small amount 2 Ad/ Int. 
Let us now determine the net shift in x due to the mo- 
tion of the wall for a walker that survives for a time r in 
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the rescaled £-plane, before it is absorbed at the right- 
hand wall (recall that t is the time in the original system, 
while r refers to the time in the rescaled £— plane, after 
the conformal mapping). 

To obtain the full solution, it is necessary to solve the 
two-dimensional diffusion equation with moving bound- 
ary conditions. Since we are mainly interested in the 
exponential tails of the winding angle distribution, we 
restrict our analysis to the limit of large times r, and 
determine the increase in x due to encounters with the 
reflecting moving wall in this limit. A Brownian walker 
that has survived for a sufficiently long time t forgets 
its initial horizontal position. The mean number of en- 
counters with the reflecting wall, and consequently the 
increase in x due to the motion of the wall, is then ex- 
pected to be simply proportional to the considered time 
interval. Assuming the validity of the central limit theo- 
rem in the limit r — > oo, the probability distribution of 
this increase Ax in x is given by 



Pa (Ax) 



1 



: exp 



(Ax - arf 
2/3V 



(19) 



The parameters a and (3 are related to the velocity of the 
wall (chirality of the defect) by a cx f3 cx v. Presumably 
Eq. (|19|) can be obtained directly from properties of ran- 
dom walks, providing the exact coefficients in the above 
proportionality. 

We can now modify Eq. (fLq) to 



Pr( x ) 



dr 



d(Ax)p/±(Ax 



3P«,i(0,t) 



di 



1 



: exp 



/2ttt 
dP a ,i{0,T 



(x - Axf 
2~7 



dr 

dr 



1 



V27tt(1+/? 2 ) 



exp 



(x — olt) 



2 1 



2r(l + /3 2 ) 



In the limit of large x it is sufficient to take the first term 
in the series expansion for P a)1 (0, r) given in Eq. (15), 
leading to 



Pr( x ) = CX P 



1 + (3 2 1 + /3 2 



(20) 



The effect of the moving wall on the winding angle dis- 
tribution is thus a systematic shift in the slopes of the 
exponential tails. For small values of chirality the slopes 
on the two sides are changed to ir/2 ± a. Due to the 
velocity dependence, these asymmetric distributions are 
clearly non-universal. At large chiralities the slopes van- 
ish as a 2 1 (3 resulting in quite wide distributions. Appar- 
ently strong chirality of a defect increases the probability 
of entanglements. 



Fig. H shows our simulation results for the winding 
angle distribution for a walk on the above mentioned 
directed triangular lattice. The asymmetry due to the 
shift is clearly visible, and the winding angle distribu- 
tion is wider than for a stationary wall. This case thus 
exemplifies the strong chirality limit discussed in the pre- 
vious paragraph. We also simulated a square lattice with 
directed bonds as indicated in Fig. ||(b) . The correspond- 
ing winding angle distribution is shown in Fig. |^. The 
distribution is again asymmetric, but not as wide as in 
the previous one, and more similar to that expected in 
the weak chirality regime. 
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FIG. 5. Winding angle distribution around the preferred 
direction for a random walk on a directed triangular lattice 
for t = 243 (dot-dashed), 729 (long dashed), 2187 (dashed), 
6561 (dotted), and 19684 (solid). The scaling variable is 
x = 2#/ln(2i). The thick solid line is the distribution given 
in Eq. §). 
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FIG. 6. Winding angle distribution for a random walk on 
a directed square lattice for t = 38 (dot-dashed), 152 (long 
dashed), 608 (dashed), 2432 (dotted), and 9728 (solid). The 
scaling variable is x = 2#/ln(2i). The thick solid line is the 
distribution Eq. (H). 
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V. WINDING ANGLES IN RANDOM MEDIA 



FLs in high-T c superconductors are pinned by oxygen 
impurities. Such pinning is quite essential to enhanc- 
ing the current carrying capacity of the superconduc- 
tor in its mixed phase |20tl . Actually, the pinning to 
the (quenched) random impurities fundamentally mod- 
ifies the properties of FLs, leading to glassy phases. The 
simplest example is again a single FL as discussed in the 
previous section. The behavior of the FL in the presence 
of point impurities can be modeled by a directed path 
on a lattice with random bond energies . In 3 or less 
dimensions, the line is always pinned at sufficiently long 
length scales. An important consequence of the pinning 
is that the path wanders away from the origin much more 
than a random walk, its transverse fluctuations scaling as 
t v , where v ps 0.62 in three dimensions, and v — 2/3 in 
two dimensions The probability of such paths 

returning to the winding center are thus greatly reduced, 
and thus the winding probability distribution is expected 
to change. 

We examined numerically the windings of a directed 
path along the diagonal of a cubic lattice (see Fig. [|(a)). 
To each bond of this lattice was assigned an energy ran- 
domly chosen between and 1 . Since the statistical prop- 
erties of the pinned path are the same at finite and zero 
temperatures, we determined the winding angle of the 
path of minimal energy by a transfer matrix method. 
For each realization of randomness this method finds 
the minimum energy of all paths terminating at different 
points, and with different winding numbers. This infor- 
mation is then updated from one time step to the next. 
From each realization we thus extract an optimal angle 
as a function of t. The probability distribution is then 
constructed by examining 2700 different realizations of 
randomness. To improve the statistics, we averaged over 
positive and negative winding angles. 

The resulting distributions are shown in figures [7] and 
||. The scaling variable in Fig. || is 20/ hit and the re- 
sults are compared to Eq. (^) that is expected in the pure 
system. A much better fit is achieved with the scaling 
variable x — 9/2^/lnt as indicated in Fig. [?]. This scaling 
form is motivated by that of self-avoiding walks, which 
in two dimensions follow a Gaussian distribution 

p flA (x = 5 A=)=-^«p (-«»). (21) 

The result of the data collapse in Fig. ^ agrees well with 
the Gaussian distribution 



p rand x = — = J — cxp (-1.5x 2 ) . (22) 
V 2 Vint/ V 7r 
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FIG. 7. Winding angle distribution for a directed path in 
a random 3-dimensional system for t — 120 (dotted), 240 
(dashed), 480 (long dashed), 960 (dot-dashed), and 1920 
(solid) . The thick solid line is the Gaussian distribution in 
Eq. @. 
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FIG. 8. Same distribution as in previous figure, but with 
x = 20/ Int. The thick dotted line is the distribution for the 
pure case given in Eq. (^). 

Directed paths in random media and self-avoiding 
walks share a number of features which make the simi- 
larity in their winding angle distribution plausible. Both 
walks meander away with an exponent larger than the 
random walk value of 1/2. (The exponent of 3/4 for 
self-avoiding walks is larger than v « 0.59 for FLs in 3 
dimensions.) As a result, the probability of returns to 
the origin is vanishingly small in the t — > oo for both 
paths, and thus the properties of the winding center are 
expected to be irrelevant. (A simple scaling argument 
suggests that the number of returns to the origin scale as 
N(t) oc l/t 1-2 ".) The conformal mapping of section || 
cannot be applied in either case: The density and size of 
impurities in a random medium becomes coordinate de- 
pendent under this mapping, as does the excluded volume 
effect. The winding angle distribution for self-avoiding 
walks in Eq. (|2~i"|) has been calculated using a more so- 
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phisticated mapping |T^,|l6|. As a similar exact solution 
is not currently available for FLs in random media, we 
resort to the scaling argument presented next. 

Let us divide the self-avoiding walk, or the directed 
path, in segments going from t/2 to t, from t/4 to f/2, 
etc., down to some cutoff length of the order of the lattice 
spacing, resulting in a total number of segments of the 
order of \nt. The statistical self-similarity of the walks 
suggests that a segment of length t/2 n can be mapped 
onto a segment of length t/2 n+1 after rescaling by a fac- 
tor of 1/2". Under this rescaling, the winding angle is 
(statistically speaking) conserved, and consequently all 
segments have the same winding angle distribution. Con- 
voluting the winding angle distributions of all segments, 
and assuming that the correlations between segments do 
not invalidate the applicability of the central limit theo- 
rem, leads to a Gaussian distribution with a width pro- 
portional to Int. This argument do es not wo rk for the 
random walks considered in section II E| IIP , since the 
finite radius of the winding center is a relevant param- 
eter. Different segments of the walk are therefore not 
statistically equivalent, as they see a winding center of 
different radius after rescaling. 

Another interesting quantitity to study is the differ- 
ence in free energy between configurations of different 
winding numbers. This provides an estimate of the free 
energy that can be gained (or lost) by a FL upon chang- 
ing its degree of entanglement. In the pure system, this 
quantity can be obtained directly from the winding an- 
gle distribution. For large t, the difference in free energy 
between configurations of winding number n and is cal- 
culated by expanding Eq. (||) for small x — Aim/ lni, and 
is given by 



F(n) - F(0) = k B T 



(lni) 2 



(23) 



(The corresponding result for absorbing boundary condi- 
tions is larger by a factor 2.) 

In the presence of quenched randomness, there ex- 
ists no obvious relation between the free energy and 
the entropy, and we therefore determined the difference 
between energies of minimal paths of different winding 
numbers. The simulation results are depicted in Fig. |^. 
Somewhat surprisingly, the difference in energy vanishes 
in the large t limit, as in Eq. (p3|). Naive arguments 
may have suggested that the energy cost associated with 
changing the winding number either saturates at a finite 
value, or possibly even grows as i 2iy_1 as in typical energy 
fluctuations JTqj . Nonetheless, Fig. ^suggests a 1/ lnt de- 
pendence on length, although we cannot rule out some 
larger power of 1/lnt. The slopes of the curves in this 
figure have ratios of approximately 1/3, 1/6, and 1/2, 
different from 1/4, 1/9, and 4/9 for an n 2 dependence. 
We coul d not collapse the data using the scaling variable 
n/VIrii as in Fig. [?|. 
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FIG. 9. Difference in minimal energy between directed 
polymers of different winding numbers, averaged over 1019 
realizations of randomness. 

Similar results are expected to hold for the windings 
and entanglements of two FLs around each other. The in- 
teractions between the FLs are small at large separations 
(very small density), and the relative distance between 
the lines behaves in the same way as the separation of 
one line from a columnar defect. In the pure system, we 
had to distinguish between repulsive, neutral, and attrac- 
tive defects. Since in the presence of point impurities the 
line does not return to the origin as often, the interaction 
with the defect is less relevant. In fact, it can be shown 
that the attraction of the columnar defect must exceed a 
finite threshold before it can pin a FL |2^] . Thus the re- 
sults of this section regarding the Gaussian distribution 
of the FL winding angle in the presence of point impu- 
rities are expected to hold even for the more realistic 
attractive columnar pins. Perhaps not surprisingly, the 
main conclusion is that the pinning to point randomness 
decreases entanglement. 



VI. DISCUSSION AND CONCLUSIONS 

Topological entanglements present strong challenges to 
our understanding of the dynamics of polymers and flux 
lines. In this paper, we examined the windings of a single 
FL around a columnar defect. By focusing on even this 
simple physical situation we were able to uncover a va- 
riety of interesting properties: The probability distribu- 
tions for the winding angles can be classified into a num- 
ber of universality classes characterized by the presence 
or absence of underlying symmetries or relevant length 
scales. 

The most "symmetric" situation is the windings of an 
ideal (Brownian) walk around a point center, described 
by the Cauchy distribution in Eq. (Q). Introduction of a 
finite core for the winding center (or a finite persistence 
length for the walk) leads to a number of exponentially 
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decaying distributions: If there is a conservation law for 
the walkers (reflecting boundaries or neutral defects), we 
obtain the distribution in Eq. (||) . Removing this conser- 
vation (absorbing boundaries or repulsive defects) leads 
to the distribution in Eq. (0) whose tails decay twice as 
fast. Constraining the end-points of the walker to the 
vicinity of the wind ing center leads to the related distri- 
butions in Sec. 



[ID 



A completely new set of distributions is obtained for 
chiral defects, where the walker is preferentially twisted 
in one direction at the winding center. These distribu- 
tions have asymmetric exponential tails, with decay con- 
stants that depend on the degree of chirality. Strong 
chirality appears to lead to quite broad distributions. A 
remaining challenge is to find the complete form of this 
probability distribution by solving the two dimensional 
diffusion equation with moving boundary conditions. 

For non-ideal walks, with a vanishing probability to 
return to the origin, the properties of the winding center 
are expected to be irrelevant. Both self- avoiding walks 
in d = 2 dimensions, and FLs pinned by point impuri- 
ties in d = 3, have wandering exponents v larger than 
1/2 and fall in this category. We present a scaling ar- 
gument (supported by numerical data) that in this the 
probability distribution has a Gaussian form in the vari- 
able 9/Vhit. Not surprisingly, wandering away from the 
center reduces entanglement. The characteristic width 
of the Gaussian form is presumably a universal constant 
that has been calculated exactly for self-avoiding walks 
in d = 2. It would be intcrsting to see if this constant 
(only estimated numerically for the impurity pinned FLs 
in d = 3) can be related to other universal properties of 
the walk. Changing the correlations of impurities (and 
hence the exponent v) may provide a way of exploring 
such dependence. 

It is likely that there are other universality classes not 
explored in this paper. The physical problem of FLs 
in superconductors provides several candidates, such as 
splayed columnar pins or a collection of pins with ran- 
dom strengths. One example, is defects with randomly 
changing chirality, which induce a change ±A9 in wind- 
ing angle each time the line returns to the defect. Such 
randomness is in fact irrelevant in the limit t — > oo, since 
it does not lead to no systematic change of the variable 
x. The induced variance in x also vanishes as 



(5Ax) 2 ~ lim \nt(^-) = 0. 

v ; t^oo I l n < J 



(24) 



The underlying problem is certainly quite rich and seems 
to call for developing some form of renormalization group 
analysis. 

It is also interesting to search for circumstances in 
which the degree of winding changes dramatically. To 
compete with the exponential tail of the winding angle 
distribution, we need an energy proportional to the wind- 
ing angle itself; a non-local quantity. A possible physical 
realization is provided by a magnetized directed polymer 



winding around a wire. A current I in the wire generates 
a magnetic field 

- x 27 X f 
W) = — . 

where c is the speed of light. Assuming that the poly- 
mer carries a uniform magnetic moment density aligned 
with its backbone, then leads to an energy 

/ dr(t') -B{r) 
'o 



proportional to the winding number n. 

In the presence of an energy e per winding, and assum- 
ing absorbing boundary conditions, the partition func- 
tion is given by 



dn exp 



en 



k B Tj 4 cosh 2 (2n 2 n/\nt) 



where n = 9/2tt is the winding number, ks the Boltz- 
mann constant, and T is the temperature. The integral 
"diverges" if t > t* = exp (47r 2 fc s T/e), a cutoff only be- 
ing given by the maximum possible winding number, pro- 
portional to t. A directed path, although free on length 
scales t < t* , is always coiled in the limit of t — > oo. How- 
ever, for a given length t, a small decrease in tempera- 
ture can induce a sharp crossover from a free to a coiled 
configuration due to the exponential dependence of the 
crossover length t* on temperature. A similarly sharp 
crossover occurs for the pinning of a FL to an attractive 
columnar defect [Q (see subsection [II D ). Both transi- 
tions are related to the number of returns of a random 
walk to the origin which scales as ln(t). We expect similar 
sharp crossovers between unentangled and braided con- 
figurations in other cases, such as several FLs winding 
around each other. 
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